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Abstract 

We show that a wide variety of generalized increasing subsequence problems admit a one parameter family 
of extensions for which we can exactly compute the mean length of the longest increasing subsequence. By 
the nature of the extension, this gives upper bounds on the mean in the unextended model, which turn out 
to be asymptotically tight for all of the models that have so far been analyzed. A heuristic analysis based 
on this fact gives not just the asymptotic mean but also the asymptotic scale factor, again agreeing with all 
known cases. 

1 Introduction 

In Prahofer and Spohn consider a certain polynuclear growth (PNG) model with stationary initial condi- 
tions, and show that it maps to the following increasing subsequence problem: 

Let t be a positive real number. Pick a random set of points in the unit square [0, 1] x [0, 1] as follows. 
On the left and bottom edges, take a Poisson process of mean t; inside the square, take a Poisson process of 
mean t^. (Thus our total mean is + 2t.) A sequence of these points is "increasing" if we have x < x' , y < y' 
whenever {x,y) and {x',y') are consecutive points in the sequence; the length of the sequence is defined to 
be the number of points. The problem is then to determine the asymptotic distribution of the length of the 
longest increasing subsequence. (Note that without the extra points on the left and bottom edge, this is just 
the standard Poisson model for increasing subsequences of random permutations 

Prahofer and Spohn then observe that the stationarity of the initial conditions can be used to show 
that the length of the longest increasing subsequence has mean exactly 2t. This fact is striking for two reasons. 
The first is that the mean in the standard model is rather complicated; it is thus surprising that a fairly 
simple extension gives rise to an exact formula for the extended mean. The second is that since adding points 
can only help the longest increasing subsequence, we conclude that 2t is an upper bound on the mean in the 
standard model. This bound is quite tight; indeed, in the standard model, the mean takes the form 2t + 0{t^^^). 
(We could also derive this upper bound from the (strictly stronger) result of that the expected length of 
the longest increasing subsequence of a permutation of length n is at most 2y^; the present method is more 
generally applicable, however.) 
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The first object of the present paper is to generalize this fact in a number of ways. It turns out, for instance, 
that the standard Poisson model admits a one-parameter family of extensions with explicit means; this further 
extends to give explicit information about the moment generating function in a neighborhood of this family. 
Also, we can replace the Poisson model by the generalized model considered in [Q (based in turn on a model of 
Johansson [pi]). 

Our other object is to explore the asymptotic relations between the extended models and the unextended 
models. It turns out that by a careful (if heuristic) analysis, we can use the moment generating function 
identities to determine not just the asymptotic mean of the longest increasing subsequence length, but also the 
asymptotic scale factor. This gives a uniform prescription for the scaling information, agreeing with the results 
of all of the cases that have so far been analyzed. 

Section 1 defines the models of interest, as well as a certain continuous limiting case. Section 2 gives a 
short, algebraic proof of the moment generating function and mean identities; this is followed by a somewhat 
more complicated, but also more enlightening combinatorial proof in Section 3. Finally, Section 4 considers the 
asymptotic consequences, giving explicit conjectures for the asymptotics of the general case. 

Acknowledgements The author would like to thank J. Baik, Y. Baryshnikov, and W. Whitt for helpful 
conversations. 



2 Extended growth models 

The model we will be considering is a generalization of a model considered by Johansson [pl| , combining the 
generalizations of and 0. We define a "parameter set" p to be a triple (i, g, r), where t is a nonnegative 

number, and q and r are sequences of nonnegative numbers with 

^g» + Xl^'<°°; (2-1) 

i i 

we will write such a parameter set as t:q/r^ omitting t \lt = Q and omitting any trailing O's from q and r. From 
the existence of the sum, we conclude that 

Q(p) := sup(gi) and R{p) := sup(ri) (2.2) 

i i 

are well-defined and attained. Given a parameter set j>, we define two functions 

H{z-p) = e*^ - zqi)-\l + zn) (2.3) 

i 

E{z-p) = e*^ + zq,){l - zn)-^- (2.4) 

i 

H{z;p) converges for \z\ < Q{p)^^, while E{z;p) converges for \z\ < R{p)^^. We will say that two parameter 
sets p+ — t^:q^ /r^ and p_ = t^:q^ /r^ are compatible if Q{p+)Q{p-) < 1, R{p+)R{p-) < 1. Note that we 
then have 

H{p+-p_) := e*^*" ne*"^«'^+'''^^e*+(«-"+'^r) _ g+q7)-i(l _ r+r;)-^{l + q+r;){l + r+qj) < ^. (2.5) 
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Let (Z+)' be a disjoint copy of Z+, and consider the set H, := [0, 1] W Z+ l+) (Z+)'. Then we associate to a 
pair p+, p_ of compatible parameter sets a random multiset C O x as follows. We let P{t) denote 

a Poisson random variable of parameter t, g{q) denote a geometric random variable with parameter g, and h{q) 
denote a random variable which is with probability and 1 with probability 

• On [0, 1] X [0, 1], we choose P{t^t^) i.i.d. uniform points, 

• On [0, 1] X i, we choose P{t^q^) i.i.d. uniform points, 

• On [0, 1] X j', we choose P{t^r^) i.i.d. uniform points, 

• On i X -we have multiplicity g{q^qj), 

• On i X j', we have multiplicity b{q^rj), 

• On i' X j', we have multiplicity g{r^r~), 

and so on, with all of the multiplicities chosen independently. 

Choose a pair of total orderings (denoted <+ and <_) on il compatible with the usual ordering on [0, 1]. A 
subsequence of a multiset M in Q x D, (a sequence of points (xi, yi) from M with no point occuring more often 
than its multiplicity) is "increasing" if we always have Xi <+ x^+i, yi <_ y^+i, subject to the further conditions 

Xi — Xi+i =^ Xi ^ Z"*" (2.6) 
y, = y,+i =^ y,iZ+ (2.7) 

In other words, the sequence must be strictly increasing along rows and columns from Z+. We then define a 
sequence Xi{M) by setting 

E ^^(^^) (2.8) 

l<i<l 

equal to the size of the longest subsequence of M which is a union of I increasing subsequences. We then define 
A(p+,p_) := A(A/(p+,p_)) (with the latter notation preferred when particular multisets are being compared). 
We recall the following results from Q : 

Theorem 2.1. For any pair p^, p_ of compatible parameter sets, A(p+,p_) is a random partition, finite with 
probability 1. The distribution of \{j)^,pJ) is independent of the choice of total orderings on VL. 

Theorem 2.2. For any pair p-^-, p_ of compatible parameter sets, we have the identity 

Pr(Ai(p+,p„) < = Hip+;p^)-' ^ueu(i) det(iJ(C/;p+)iJ(C/t;p_)). (2.9) 

Remark. As in [Q, this is a formal integral, defined by analytic continuation from the region Q{p±) < 1. 

For our purposes, we need to extend the model slightly, by adding a special row and column to the random 
multiset. Extend to il^ by adding a new element, denoted S, and extend the total orderings so that S is 
smallest in both orderings. Then we define a new random multiset Af (p+,p_; a+, a_) as follows. 
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On ft X n, we take M{p+,p^), 

On [0, 1] X {S}, we choose P(a_t+) i.i.d. uniform points, 
On we have multipUcity 



• On («', S), we have muhipHcity g(a_r^), 

and so on, and allow increasing subsequences to be weakly increasing in the new row and column. Note that 
the point (S, E) has multiplicity fixed at 0; otherwise, the new model would simply be a special case of the old 
model. By the argument in we have: 

Theorem 2.3. For any compatible pair p+, p_ of parameter sets and any a+, «_ with a+R{p-) < 1 and 
a-R{p-^) <1, we have 

Pr(Ai(p+,p_; a+,Q;_) < = 

E[a+]p_Y^E{a_]p+y^H{p+]p_Y^[Di{p+,p_;a+,a_) - a+a_ «+, "-)], (2.10) 

where 

Di{p+,p^-a+,a^) - Eyec/(0 det((l + a+U){l + a^U^)H{U-p+)H{U^ ;p^)). (2.11) 

Aside from intrinsic interest (|7j), this new model in principle can give us some information about the original 
model Ai(p+,p_), since we have the coupling 

Ai(M(p+,p_)) < Ai(M(p+,p_;a+,a_)). (2.12) 

Thus it is particularly interesting that, as we shall see, there is an exact formula for 

E(Ai(p+,p_;a,a-i)) (2.13) 

whenever R{p+) <a < R(p^)^^. 
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We will also consider a continuous limiting case of the above model, combining the exponential limit of 
and the heavy-traffic limit of queueing theory Q (studied in the present context in ||l^ and ^). 

The parameters for the continuous model consist of a pair of sequences p^eRUjoo}, a nonnegative real 
u, and two numbers a± G M U {oo}, subject to the convergence constraint 

^ < oo, for z < inf (p±), (2.14) 

and the compatibility constraints 

inf(p+) + inf (p-) > 0, inf(p+) + a_ > 0, a+ + inf (p") > 0. (2.15) 
We omit any trailing go's from p*. Define an infinite matrix Ay, < i, j hy 

fo (*,j) = (0,0) ^ ^ 

Ai, = < (2.16) 
(z,j)^(0,0). 
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with the convention that = a±, and let M be a random matrix fiUed with independent exponential random 
variables, such that Mij has mean A^ . Also, for < i such that < oo, let Bi be a Brownian motion on [0, 1] 
with 

E{B,{x)) ^ -up+x Y&r{Bi(x)) = ux; (2.17) 

by convention, — a±. Then we define a random sequence a+, a_) in terms of increasing paths 

in N X (N 1+) [0, 1]), where the contribution of an interval i x [xi,yi] is Bi{yi) — Bi{xi) and the contribution of a 
point i X j is Mij. When w > 0, we require that the union of k increasing paths used to define Xk have total 
Lebesgue measure fc in N x [0, 1], and that the paths contain no point in i x [0, 1] when pf — oo. If no such 
path exists, we set Xk = — oo. 

Theorem 2.4. For any valid choice of p^ , u:p^ , a+, a_, of compatible parameter sets, x(/9^, a+, a_) 

is a random nonincreasing sequence m MU {— oo}, and is nonnegative if u — 0. The distribution is invariant 
under reordering of the sequences p^ . 

Proof. That the sequence x is nonincreasing, invariant under reordering, and nonnegative when u = follows 
from the fact that it is a scaled limit of A (see below); it suffices therefore to show that xi < oo with probability 
1. We first observe that 

Xi(p+,u:p";a+,a_) < xi(p+, uip"; a+, oo) + xi(p+, oo, a_). (2.18) 

Since decreasing a± cannot decrease xij we conclude that it suffices to prove finiteness when — inf(p^) < a± < 



inf(/3^) and a+ + a_ > 0. But then Corollary 3.4 below applies, expressing 'E(e''°-+^°'-^^^ ) as an infinite product. 
Moreover, the conditions on p^ suffice to force convergence of this product, and thus E(e'^"++°"^'^i ) < oo. But 
this immediately implies xi < oo, as desired. □ 

As we alluded to above, this is a limiting case of the discrete model: 

_ X{/e-''^/\ut^:/e-'>-/';e-^/\e---/*) - ut^ 
X[P ,u:p ;a+,a_)= hm , (2.19) 

with the limit taken in the sense of distribution. This corresponds to the facts that if x is an exponential 
random variable of mean 1/m, then [jj is a geometric random variable of parameter e~™/', and that as the 
parameter tends to infinity, Poisson processes converge (with proper scaling) to Brownian motion. 

The main significance of the continuous model is that it contains two classical matrix ensembles as special 
cases. Let 0„ denote the finite sequence consisting of n copies of 0. We obtain the Gaussian Unitary Ensemble 
from x(On, oo, oo) that is, x(On, 1:; oo, oo) is distributed as the (ordered) eigenvalues of an n x n Hermitian 
Gaussian matrix (extended to an infinite sequence by adding — oo). (This is essentially proved in [ p^ ; to be 
precise, they show that one obtains traceless GUE under a constraint equivalent to J^i ^i(l) ~ 0; it follows easily 
that without this constraint, one obtains ordinary GUE.) Similarly, we obtain the Laguerre Unitary Ensemble 
from x(On+ , 0„_ ; oo, oo); that is, the distribution of the singular values of a x n_ complex Gaussian matrix. 
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3 An algebraic proof 

Let p_ be a pair of compatible parameter sets. 

Lemma 3.1. Di{a^, Z?/ (p+ , p_ ; a+ , q;_ ) is a polynomial, satisfying the identity 

Di{a+,a-) = a;'). (3.1) 

Proof. This follows immediately from the corresponding fact for det((l + «+[/)(! + a^W)). □ 

For a+ < R{p-)^^, a_ < define L(a+,a_) = Ai(p+,p_; «+, then 

Theorem 3.2. For R{p+) < a+ < R{p-y^ and R{p-) < a- < R{p+y^ , 

E((a+a_)-^("+^"-)) - Eia+;p^)-'Eia^;p+)-'Eia^';p+)EiaZ';p-) (3.2) 



Proof. By Theorem 2.3 above, 

Pr(L(a+,a_) < I) ^ E{a+;p^)-^E{a-;p+)-^H{p+;p^)-^[Di{a+,a^) - a+a^Di^i{a+,a^)] (3.3) 

= {a+a^yEia+;p^)-'Eia-;p+)-'H{p+;p^)-'[DiiaZ\a^')-Di^i{aZ\a-% (3.4) 

Then 

^ (a+a_)-'Pr(L(a+,a_) = 

{a+a-)-^[Pr{L{a+,a^) < I) - Pr(L{a+, a-) < I - 1)] (3.5) 

0<i<fc 

(a+a_)-'=-i Pr(i(a+, a_) < fc) + "+"~ ~ ^ V {a+a^y^ Pr(L(a+, a_) < /) (3.6) 

^ 0<l<k 

Eia+;p^)-'E{a_;p+)-'H{p+;p_)-' 

(3.7) 



0<i<A: 



= £;(a+;p_)-ii;(a_;p+)-ii/(p+;p_)-i[i?fe(a:\a;i)-(a+a_)-ii?fe_i(a:\a;i)] (3.8) 

= E{a+;p^)-'E{a^;p+)-'Eiay;p+)E{aZ';p-) Pr{L{aZ\a^') < k), (3.9) 

where the last step is valid since aZ^ < R{p-)^^ and < R{p+)^^. 

The theorem then follows by taking the limit k —^ oo. □ 

Taking a limit as ^ 1, we obtain: 

Corollary 3.3. Whenever R{p+) < a < R{p^)~^ , 



„i aE {a;p^) a ^E {a ^■,p+) 
E i(a, a ) = — + — -— — . 

E{a;p-) E[a ^\p+) 



(3.10) 
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Remark, li p = t:q/r, then 



E{a;p) ^ 1 



We also observe that 



aE'(a;p-) a''^ E' (a''^ ; p+) d , E(a;p^) 
' — a— log ■ 



(3.11) 



(3.12) 



E{a;p-) E{a ^;p+) da E{a ^;p+) 

In the continuous limit, we write X(a+,a_) := xi(p^i ^^^P ; a+j o-)- Taking the appropriate limit gives: 

Corollary 3.4. Let p+ and u:p^ be compatible continuous parameter sets. Then whenever — inf(p^) < a+ < 
inf(p+) and — inf(/9+) < a_ < inf(/3^), 



g(-g(a++a_)X(a+,a_)-j ^ ^u{al-al)/2 "Q 



n 



Whenever — inf(p ) < a < inf (p+), 

E{X{a, -a)) = -ua + J2 + J2 

4 A combinatorial proof 

Fix parameters as in the previous section, and set 

N+ = |M(p+,p_;a+,a_)n ({S} x n)\ 



iV_ = \M{p+,p,-a+,a-)n{n X {S})| 



Then we observe 



and for a+ = a_ 



E((a+a_)-^^+) 



E(iV+) 



i?(aI^P-) 
£;(«+; p_) 

E'{a-,p-) 
' ^(a;p-) 



E((a+Q:_ 



E(Ar_ 



£:(a_;p+) ' 

E'{a-^\P+) 
E{a-^;p+) • 



So we can restate Theorem and Corollary 3.3 as 

E((a+a_)-^("+'"-' 



(3.13) 



(3.14) 



(4.1) 
(4.2) 

(4.3) 
(4.4) 



(4.5) 
(4.6) 



E((«+a_)-^+-^-) 
E(L(a+, a_)) = E(iV+ + 7V_) 

We give a direct proof of this fact, for a+a- < 1: 

Proof. Let a' < Then a+a' < 1, so we can extend M{p^,p^] a+, a') by adjoining (E, S) with multiplicity 
Nq of distribution ^(a+Q;'); denote the resulting random multiset by M' . But then by Theorem 2A, we can 
change the total ordering <_ so that S becomes maximal instead of minimal. We then find that 



S X 17, (S,I]),rj X s 



(4.7) 
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induces an increasing subsequence of M' with respect to this new ordering; thus 

No + N+ + N_ <\i{M'). (4.8) 

On the other hand, this is the only maximal increasing subsequence that passes through (S,5]); any other 
maximal increasing subsequence can have size at most Nj^ + N- + Ai(M(p+,p_)). We thus find 

Ai(M') < Ar+ + AT. +max(Ai(M(p+,p_)),7Vo). (4.9) 

But Ai(M') is distributed as A^o + •^i(Af(p+jP-; 0+, a')) (since before the reordering every maximal increas- 
ing subsequence passes through (S,S), and Nq is independent of Ai(M(p+,p_; a+, a')). So if we take the 
expectations and subtract/divide by the contribution of A^O) we find that we need only show 

lim E(max(Ai(p+,p_),iVo))-E(7Vo) = (4.10) 

a'— »(l/c«)~ 

lim E((a+a_)-^°)-iE((a+a_)-'"'«(^^(f+'P-)'^'')) = l. (4.11) 

Lemma 4.1. Let X he a nonnegaMve-integer-valued random, variable with finite first moment, and let Y be an 
independent geometric random variable of parameter t. Then 

lim E(max(X, Y)) - E(F) = 0. (4.12) 

Similarly, for s < 1, ifE{s~^) is finite, then 

lim E(s-^)-i E(s- ^'^^X'Y)) = 1. (4.13) 



Proof. 



Similarly, 



and thus 



E(max(X, Y)-Y)= E((X - Y) Pr(F < X)) (4.14) 
= E(X) + ^(E(t^)-l) (4.15) 
^ 0. (4.16) 



j,^g-max(x,y)) = E(s-^ Pr(F < X) + Pt{Y > X) E(s-'^|y > X)) (4.17) 
= E(s-^) + fc-^E((Vs)^). (4.18) 



lim E(s-'^)-i E(s-'"^(^''^)) = lim j \ E(s-^) + E((i/s)^) (4.19) 

f— s-s >5 Syt — Ij s{t — 1 j 

= 1 (4.20) 

□ 
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The theorem then follows from the following lemma, since 



a+a- > R{p+)R{p-) (4.21) 
by assumption. □ 

Lemma 4.2. For all z > R{p-^-)R{p^), Ei{z^-^^'-p+'P~^) is finite. In particular, since 1 > R{p+)R{p-), 
Ai(p+,p_) has moments of all orders. 

Proof. An increasing subsequence in can pass through a point on a strict row or column at most 

once; thus Ai(M(p+,p_)) is unchanged if we remove any excess multiplicity in those rows and columns. Let 
M° be the resulting multiset, then 

Ai(M(p+,p_)) = Ai(Ar) < |M°|. (4.22) 

It will thus suffice to prove that E(z^l*^ ') < oo. But the moment generating function of \M°\ is /(z)//(l), 
where 

f{z) = e-*^*"- []e*^9r-e*+'T^e«^'^e'^^"^ ^(1 - r+rjz)-\l + g+r7z)(l + r+qjz){{l - q+qj) + q+ qj z) 

(4.23) 

This product converges to an analytic function with no pole inside the open disc \z\ < (i?(p_|_)_R(j<_))^^, and 
thus the result follows. □ 

5 Asymptotic consequences 

Since Ai(p-|-,p_; «+, q;_) is nondecreasing in Q!_|_ and a-, we obtain the following bound: 
Theorem 5.1. For any compatible parameters p 4^, p-, 

E(Ai(p+,p_)) < inf E(Ai(p+,p_;a,a-i)). (5.1) 

R.{p+)<a<R(p-)-^ 

For instance, in the purely Poisson case, p+ = p_ — t:/ , we find 

E(Ai(p+,p_)) < inf (a + a"i)t = 2i. (5.2) 
This bound is remarkably tight; indeed, we have |^ 

Ai(p+,p_) = 2t-0(ti/3+.). (5.3) 

This suggests the following conjecture: 

Conjecture 5.2. Fix parameters p+ , p-, and define 

m{a;p+,p-) = E(Ai(p+,p_; a, a"^)). (5.4) 



Then 

lim Xi{p'l,p"i) — inf m{a;p+,p-), (5-5) 

n->oo _R(p+)<Q<fl(p_)-i 

with probability 1, where 

{t:q/r)" := (nt):q"/r" (5.6) 
g" := . ..qi,q2,q2, . . .^2, • • • • (5.7) 

Remark 1. Roughly speaking, this is an analogue of the law of large numbers. As such, it can most likely 
be strengthened considerably (considering different sequences of parameter sets than just p", p"). See, for 
instance, the result of [[l6[. 



Remark 2. The existence of the limit (5.5) follows from superadditivity and the bound ( ^.1| ). 

This has been verified in a number of special cases (see below). In each case, we in fact find that 

Ai(p",p") - 



(5.; 



ni/3 

converges to a limit distribution. 

Fix parameters p+, An increasing subsequence of Af (p+,p_; a, a~^) cannot include points from both 
{S} X Vl and Q. x {S}. We would thus expect that for a "large" , the typical longest increasing subsequence will 
avoid r2 X {E} entirely. In particular, we would expect 

E(Ai(p+,p_;a,a-i)) ~ E(Ai(p+,p_; a, 0)) (5.9) 

whenever A^_|_ 3> N^. For asymptotic purposes, this condition is simply that a > a+, where d+ minimizes 
m(a;p+,p_). (We also define a_ — (a+)^^, which of course minimizes m(a;p_,p+)). 

In particular, this tells us that d± is a critical point; if a+ < q:_|_ and Q!_ < (5_, we have 

E(Ai(p+,p_;a+,a_)) ^ E(Ai(p+,p_; a+, a_)), (5.10) 

while if either is greater, the mean is determined by the dominant parameter. 

This behaviour is, of course, confirmed by the analysis of Q, in which the asymptotics for general a± are 
determined for the Poisson case p = p' = t:/ and the Johansson case p = p' = j ^fq^ (where y^" is the finite 
sequence consisting of n copies of ,Jq). In both cases, we obtain the same behaviour near the critical point. 
This suggests that for general parameters there should exist constants /i, cr, and a± so that the following holds: 

If we fix vjj^ , W- , and define 

a±=a±exp( ^^^), (5.11) 



then as n — *■ cxD, 



converges to the distribution H{w+,W-) (0). 

We recall the following information about the distribution H(w^^,w 



Ai(p",p";Q+,a-) - fJ-n ^^^2^ 
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Lemma 5.3. Let X be distributed as i7(w+,w_). Then E(exp(2(w+ + w^)X)) = cxp(|(wi|. + w't))- If 
w+ = — w_ = w, then E(X) = iw^ . 

Only the latter equation was actually shown in |Q, but essentially the same calculation gives the first 
equation as well. Furthermore, in the cases that have been fully analyzed, this is precisely the analogue for H 
of Theorem 3.2 and Corollary 3.3 above. This suggests that we compare the results in general. 

On the one hand we compute 

log(E((a+a_)~^i(f^'f-;"+'"-))) = 2/in2/='(^ + — ) + log(E(exp(2w+Xa/a+) exp(2w_Xa/a_))), (5.13) 

(7+ (J- 

with X distributed according to H{w+, w_) in the limit. Thus to retain the analogy, we must have (7+ = (t_ = a. 
On the other hand, we have: 

log(E((a+a_)-^i(P?'f-'"+'"-))) = 

2nV3(,,)(5^)!!!±±^ + 2nV3(,2,)(5+)!^:l_^ + l(^3^)(«+)!^^i±^ + 0(^-^3). (5.14) 

(T (T^ 3 O"^ 

where {9f)[z) — z-j^f{z) and g{z) = logE{z;p^) ~~ log E[z^^ ; p-^-) . Thus, in particular, [6g)[z) — m{z:p^,p). 
Comparing the asymptotics, we find 

a = ((^3^)(a+)/2)i/3 ^ (ci>"(5+;p+,p_)/2)i/3. (5.15) 

Similar considerations (based on part (iv) below) give us the scale factors for a+ > a^^-, thus giving us the 
following conjecture (see for the definitions of the limiting distributions): 

Conjecture 5.4. Fix parameters p+, define a± as above, and further define 

fi = m{a+;p+,p^) ct = (a+m"(a+;p+,p_)/2)^/^ (5.16) 

fi+(z) = m{z\p+,p-) a+{z) = {zm'{z;p+,p-)y^^, (5.17) 

fi^{z) ^ m{z;p-,p+) a-{z) ^ (zm'{z;p^,p+)y/'^. (5.18) 

Assume that a± ^ {0, 00}, a > 0, cr+(Q:+)^ > for Q!_|_ < q;+ < R{p^)^^ , and (T_(a_)^ > for a_ < a_ < 
R{pj^)^^ . Then we have the following limiting distributions as n — )■ 00; 

• (i) If < Q!+ < ci_|- and < q;_ < (5„ are fixed, then 

Ai(p",p";a+,a-) - fin 

Near the critical point, set w± by a± — a±e'xp{—2w±/ an^^^). 

• (ii) If w± and < aip < cizp are fixed, 

Xi{p+,pl;a+,a-) - fin 



Fgve (5.19) 



1/3 ^ G{w±). (5.20) 
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• (Hi) If w+ and are fixed, 



H{w+,w-). (5.21) 



Finally (the Gaussian regime), let a^. and a*L be fixed such that a!]. > d± and /z+(a*]_) — /i_(a'L), and set 

• (iv) If x± and < < a^j- are fixed, 

^ ^ > N{x±,(j±{a^) ). (5.22) 



• (v) If and x^ are fixed, 

„l/2 



mayi{N{x+,a+ialf),N{x-,a-ia°_f)). (5.23) 



Remark 1. This, in turn, is analogous to the central limit theorem, so again can probably be strengthened 
considerably (although not nearly to the same extent as Conjecture most likely can). In particular, it is 
presumably sufficient for the parameters a±, w±, x± to tend to limits as appropriate, rather than simply be 
fixed. 

Remark 2. We recall that 

/i+(a) =E(iV+) + E(iV_) (5.24) 

at a+ = a, q;_ = a~^. Similarly, 

a+{a) = Var(iV+) - Var(7V_) (5.25) 

Remark 3. The analogous conjectures for models of the other symmetry types (Q, ||^, [||) are straightforward. 
We note in particular that when = p- = p, we find am'{a) — —a^^m'(a^^), and am'(a) > whenever 
a > 1. So the hypotheses of the above conjectures hold in such cases, with a± — 1. 

In the continuous limit, we make a similar conjecture; the main difference is that the model is nonincreasing 
in the parameters, not nondecreasing, so the -Fgue region is now a± > d±. The scale factors are: 

fi^m^{d+;p+,u:p-) a ^ im'^{a+; p+ ,u:p-)/2y/^ (5.26) 

p+{z) = mc{z; p'^ ,u:p^) cr+(z)^ = -to^(z; p+, mp"), (5.27) 

P-{z) = p+{-z) a-{zf = -a+{-zf (5.28) 



where we define 



mc{z; p'^ ,u:p ) = — tiz + ^ 1" ^ (5.29) 
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Near the critical point, we take = + 2w^ /an^^^ , while in the Gaussian regime, we take = — 
x±/CT±(a±)ni/2. 

As remarked above. Conjecture 5.4 was proved in |Q (with the exception of parts (iv) and (v), which are 
straightforward using the argument in section 7 of for the cases p± = t:/ and p± ~ /v^"- The only other 
known results are for the case a± = 0; the references in the following examples refer to this case alone. 

Example 1. If we take p± = t:/, we have 

fi±{z) = {z + z-^)t, cr±(z) = ((z - z"i)i)i/^ a=l, p = 2t, a^t^^^. (5.30) 
This corresponds to the classical case of increasing subsequences of random permutations, studied in |^ . 
Example 2. If we take p± = /(-v/?)"*, with n^/n^ tending to a constant, we have 



lJ-+{z) 



1-z^ z-^' 



a+izY 



and thus 



(1 - z^y [z 



(5.31) 



/qnZ+ JtH^ 



1-9 



1-9 



This model was analyzed in |[ll[| , along with the continuous (Laguerre) limit, in which case we have 



M+(2) 



1 r Tij 



2'l + 2z 1-2Z-' 



a+{zf 



4n I 



4n_ 



and thus 



Example 3. If we take p+ = /I", p_ = t:/ , we have 

pL+{z) = Zt + 

and thus 



(l + 2z)2 (l-2z)2 



z- 1' 



cr_|_(z)^ = zt - 



zn 



(z-iy 



a+ = l + VW*, p- = t + 2Vtn, a ^ Vtn^'^^^{l + \/^) 



2/3 



(5.32) 



(5.33) 



(5.34) 



(5.35) 



(5.36) 



This model, corresponding to weakly increasing subsequences of random words, was studied in |10[ and |lq| . In 
the continuous limit (corresponding to the n x n GUE p8|), we have: 



/X = 2-\/n, a 



,-1/6 



(5.37) 



These are precisely the scale factors required to make the largest eigenvalue of an n x n Gaussian Hermitian 
matrix tend to the limit -Fque- 
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Example 4- If we take p+ = (v^)"+/, = (^Z?)" /> with n-f-/n^ tending to a constant, we have 

= -r- = H = = -; —-r - —-r. (5.38) 

Here we have three cases. If qn^ > n_, then q:_|_ = 0, and if qn^ > ri+, then 5+ = oo; in either case, 
the above conjectures do not apply (indeed, in those cases one expects the hmiting distribution to be atomic, 
Ai(p+,p_) = min(n_|_, Otherwise, 



o^+ = —== '-^=, lJ- = — , cr = 



n+ - y/qnZ 1-q 1-q 

(5.39) 

The mean in this model was derived in p5|; the refined asymptotics of a symmetrized version was studied in 



Example 5. If we take p+ = (y^)"+/, P- — /(v^)"^i with n-^-/n^ tending to a constant, we have 

There are two cases. If qn^ > n_, then — 0, and the above remark applies. Otherwise 

— -y/W+ 2^qn+n^+q{n^ - n+) (.^^H'^-j '(1 - ^ — ) '(1 + ^ — J ' 



(5+ = — ^ • /i = — ^ . cr = 

y^rTjT + ^gn_ 1 + 9 ' 1 + 9 

(5.41) 

The mean in this model was first derived in p7| ; the fluctuations have been analyzed in section 5 of ]lOt . 
Example 6. If we take p+ = 1"/, p- = t:/ , we have 

Th ZTl 

^i+{z) = zt+—-, <T+{zf = zt - -——^ (5.42) 
z + 1 (z + 1)^ 

If n < t, then (5_|_ = 0; otherwise we have 

(5+ = y^-l, ^ = 2Vt^-t, G = (tnfl^{\- ^ftfTifl'^. (5.43) 

This corresponds to strictly increasing subsequences of random words, and was studied to a small extent in [p^ . 

Note that the pathological case a+ = (resp. (5+ = oo) can only occur if all the rows (resp. columns) are 
strict. 
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